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Abstract. We have examined data from pulsars B0950+08 and B0329 +54 for evide nce of temporally coherent 
radiation using the modified coherence function (MCF) technique of lJenet et alJ (1200 We consider the influence 
of both instrumental bandpass and interstellar propagation effects. Even after removal of the effects due to the 
instrumental bandpass, we detect a signature in the MCF of our PSR B0329+54 data which is consistent with the 
definition of a coherent signal. However, we model the effects due to interstellar scintillation for this pulsar and 
show that it reproduces the observed signature. In particular, the temporal coherence time is close to the reciprocal 
of the decorrelation bandw idth due to d iffractivc scintillation. Furthermore, comparison of the coherence times of 
three pulsars reported by IJenet et alJ (l200lT) with their expected diffractive decorrelation bandwidths suggests 
that the detection of coherence in these pulsars is also likely a result of interstellar scintillation, and is not intrinsic 
to the pulsars. 
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1. Introduction 

Only a f ew months after the discovery of the first pulsars 
in 1967 l|Hewish et alJll968l) it was realise d they were in 
fact neutron stars fcoldlll968l[Pacinilll968lb However, the 
emission process, responsible for the radio-frequency radi- 
ation in pulsars, remains uncertain. One characteristic of 
the pulsar radiation is a very high brightness temperatur e 
which can reach 10 30 K (e.g. iManchester fc Tavlorlll977j) . 
which implies that it is generated by a coherent mecha- 
nism. It is therefore important to identify charact e ristics 
of coherent emission in pulsar radiation. iRicked lll975ft 
suggested that the observed radiation is well described 
by amplitude modulated Gaussian noise, which contains 
no i nformation ab out the coherent nature of the radia- 
tion. [Corde3Jl97(|. later generalized this to non-Gaussian 
shot noise. Irlankins et al.l (|2003il have seen individual shot 
nano-pulses in the giant pulses from the Crab pulsar, sug- 
gesting that the intrinsic structure is less than 2 ns. 

IJenet et alJ l)200l|) (hereafter JAP) claim to have de- 
tected the existence of coherent non-Gaussian radiation 
on 100 ns time scales for pulsars B0823+26, B0950+08 
and B1133+16 in observations made at the Arecibo ob- 
servatory. They used their "modified coherence function" 
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(MCF), defined in Eq. below, to find statistics inconsis- 
tent with amplitude-modulated Gaussian noise in the volt- 
age time series from observations of these three pulsars. 
They show that a coherent model for pulsar radiation can 
account for the observed statistics. However, if the MCF 
of a time series of pulsar radiation shows statistics of a 
non-Gaussian nature, it does not yet prove that the pulsar 
radiation itself contains non-Gaussian statistics. We have 
to consider the effect of the ISM on pulsar radiation and 
the instrumental effects on the signal after detection. JAP 
assert that scintillation does not influence the statistics of 
the signal measured by the MCF. In particular, they state 
that if the statistics of the intensity fluctuations are well 
described by Gaussian noise, the MCF ought not to ex- 
ceed zero as a result of propagation effects. Scintillation, 
however, can give rise to a quasi-periodic fluctuation in the 
frequency power spectrum (see Fig. [2b) . Such fluctuations 
in the frequency domain can be expected to influence the 
autocorrelations from which the MCF is constructed. 

We consider here an independent data set for pulsars 
B0950+08 and B0329+54 to try to confirm the detection 
of coherence. In Sect. [21 we describe our observations and 
the construction of the MCF. In Sect. [31 we present an 
analysis of our data and investigate whether scintillation 
and instrumental effects do influence the MCF. Here, we 
also present a numerical experiment that shows that scin- 
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Table 1. Pulsar Characteristics and Analysis Parameters 



Pulsar 


B0950+08 


B0329+54 


Period 


0.2531s 


0.7145 s 


Period derivative 


8.6 x 10~ 17 


2.04959 x 10~ 15 


Dispersion measure 


2.9702 pc cm" 3 


26.7765 pc cm -3 


Centre frequency 


382 MHz 


328 MHz 


Bandwidth 


10 MHz 


5 MHz 


Number of Pulses 


2 320 


5 200 


Size of window 


102.4 /xs 


204.8 fis 


Number of windows 


16 


16 


Time-resolution 


50 ns 


100 ns 



filiation quantitatively explains the coherent features ob- 
served by JAP and compare their values of the coherence 
time with the inverse of the diffractive decorrelation band- 
width for each pulsar. In Sect.^we discuss our results and 
present our conclusions. 

2. Observations and analysis 

For our analysis, we use data from pulsars B0329+54 and 
B0950+08. These are bright pulsars and it is therefore 
possible to detect single pulses with the high signal-to- 
noise ratio required to find a clear signature of coherent 
radiation. 

These data were taken using the Westerbork Synthesis 
Radio Telescopes ( WSRT) with its pulsar backend, 
PuMa IjVoute et alJ 12002'). In its tied array mode the 
WSRT is equivalent to a single dish with a diameter of 
94m and has a gain of 1.2K/Jy. PSR B0950+08 was ob- 
served on 28 April 1999 at a centre frequency of 382 MHz, 
with a bandwidth of 10 MHz. PSR B0329+54 was ob- 
served on 13 August and 11 September 1999 at a cen- 
tre frequency of 328 MHz, with a bandwidth of 5 MHz. In 
both observations a 10 MHz band was Nyquist-sampled 
in 2 linear polarisation channels. For PSR B0329+54 the 
5 MHz band was formed by digitally filtering the data us- 
ing a finite impulse response filter. After sampling, the 
data were 2 bit digitised. The high time resolution of 50 ns 
and 100 ns for pulsars B0950+08 and B0329+54, respec- 
tively, is necessary to study radiation which is expected to 
contain featu res with a coheren ce time of a few hundred 
nanoseconds (jjenet et al.ll200lj) . 

We remove the effects of i nterstellar disp ersion by 
means of coherent dedispersion l)Hankinslll97llh Finally, 
for both pulsars, we take two sets of small successive time 
windows for each individual pulse. One set contains small 
windows centered on the peak of the average pulse pro- 
file, the other set contains a region outside the pulse, 
which is used for system plus sky noise corrections. The 
size, number of windows, and other parameters can be 
found in Tabled To make sure that the periodicity of mi- 
crostructure does not show up in our results as coherence, 
we chose a window smaller than the typical microstruc- 
ture time-s cale for both pulsars , which is 0.17 ms for PSR 
B0950+08 faickett et al.lll975l) and 0.6 - 1.5 ms for PSR 
B0329+54 llLange et alJ IlE 



From these windows we calculate the MCF for real 
voltages. The MCF tests the fourth moment of the signal 
against the square of the second moment and is given by 



M s (A</>) 



_ (Cis(M)) 

(Ci s (0)) 



(Cs(0)) 



(1) 



where S is the time series of real voltages, Is — S ■ S is 
the intensity, while Cs(A<^>) and C/ s (A0) are the auto- 
correlations of S and Is, respectively. The angular brack- 
ets denote averaging over different windows and different 
pulses. In the analysis of JAP, S is a time series of com- 
plex voltages. Although the MCF for complex voltages is 
different from the MCF for real voltages, physically, they 
are the same. 

The MCF is sensitive to phase relationships between 
measurement points separated in time, called temporal 
coherence. The two polarisations are also separated and 
treated as different pulses. We thus obtain 32 time series 
from each pulse. 

2.1. Calculating the MCF 

The MCF is sensitive to coherent features present in any 
time series. For pulsar radiation, such features could re- 
sult when the radiation is emitted in bunches of coherentl y 
radiating particles (e.g. iRuderman fc Sutherland! Il975]h 
However, if there is extensive incoherent addition of these 
bunches, t he resulting ra diation will contain only Gaussian 
statistics l|CordesllT97f^ . Because of the way the MCF 
is constructed (see Eq. (QJ), it has a value of zero for 
Gaussian noise (if applied to complex data) and falls be- 
low zero with increasing delay for amplitude modulated 
Gaussian noise. Therefore, if the MCF of a time series 
becomes larger than zero, it contains statistics of a non- 
Gaussian nature, which might give information on bunch- 
ing of coherently radiating particles. To find such statis- 
tics, we wish to calculate the MCF for a noise-subtracted 
pulsar signal. 

In order to calculate the second term of the MCF we 
note that the recorded on-pulse voltage time series V(t) 
contains both the pulsar signal and system plus sky noise. 
It can therefore be expressed as 



V(t)=S(t)+N(t), 



(2) 



where S(t) is the pulsar signal and N(t) is the system plus 
sky noise, which is contained in the off-pulse signal. Using 
the voltage time series of the on- and off-pulse, the av- 
erage autocorrelation (AAC) of the pulsar signal and the 
AAC of the intensity of the pulsar signal are calculated as 
follows. First, one window containing one polarisation of 
the on-pulse of the pulsar profile is Fourier transformed 
to the frequency domain. The value at each frequency is 
multiplied by its complex conjugate and inverse Fourier 
transformed back into the time domain (correlation theo- 
rem) . This results in the temporal autocorrelation of one 
window of on-pulse signal. The autocorrelations of all the 
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(a) 



(b) 



Fig. 1. Our results for the MCF as a function of phase delay calculated for pulsars B0950+08 (a) and B0329+54 (b) 
before bandpass correction, obtained from 2 320 and 5 200 time series of real data obtained with WSRT, respectively. 
The time resolution is 50 ns for PSR B0950+08 and 100 ns for PSR B0329+54. Note that the peak at zero phase delay 
has been set to zero. 



windows in the total time series of on-pulsc signal (con- 
taining both polarisations) are then averaged together to 
give the AAC of the on-pulse, (Cy (Acj))). The same calcu- 
lation is made for the off-pulse signal to give (C,/v(A</>)}. 
The AAC of the pulsar signal without the system and sky 
noise, is given by 



(Cs(A^)} = (C v (A<f>)) - <CW(A0)). 



(3) 



This AAC is divided by the average pulsar intensity of the 
total time series, which is simply the average intensity of 
the on-pulse minus the average intensity of the off-pulse, 



(Is) = (Iv) - (In). 



(4) 



Squaring the result, we obtain the second term of the 
MCF, (C S (A0)) 2 /(C S (O)) 2 , where (C s (0)) = (I s ). 

To obtain the first term, we calculate the AAC's of the 
intensity of the voltage time series of on-pulse signal and 
off-pulse signal. The AAC of the intensity of the noise- 
subtracted pulsar signal is given by 

(C Is (Act>)) = (C Iv (&ct>))-(CiA&<t>)) 

-2[(CV(0))-(CW(Q))]-(C!y(0)> 

-4 [(CV(A£)) - (CV(A$>] • (CV(A0)). (5) 

The first term in the MCF results by dividing this AAC 
by the average of the square of the intensity of the total 
time series, which is given by 

(C 7s (0)) = (C Jv (0))-(C Jw (0)) 

-6[(C V (0))-(C N (0))]-(C N (0)). (6) 

The MCF is then calculated according to Eq. . 
3. Results 

Our MCF for PSR B0950+08 is shown in Fig. HJl. The 

MCF starts at a value of zero and falls with increasing 



A0, which is due to the non-zero slope of the average 
pulse profile of the pulsar. Our MCF for PSR B0329+54 
is shown in Fig. ^p. It shows that there is a broad ex- 
cess reaching a Acj) of 35 ^s. There is also a large peak 
extending over the second and third bins. 

The peak in PSR 0329+54 is at a Acj) that is the re- 
ciprocal of the bandwidth. This makes us suspect that 
the frequency-dependent gain of the system, the bandpass, 
can influence the MCF. We can understand the presence 
of a peak in the MCF due to the bandpass as follows 1 . 
Define S(v) as the Fourier transform of the voltage signal 
S(t) and define the bandpass as S(v)S* (u). The bandpass 
contains modulation with a width equal to the width of 
the bandpass, Ai/b p and also modulation with a width of 
about one third the width of the bandpass (see Figs. Et 
and d). This implies that S(t) and S(t) 2 contain a peak 
with a width on the order of 3/A^b P - The autocorrela- 
tion of the signal Cs(At) and the autocorrelation of the 
intensity Cj s (At) will then contain a peak at At — and 
decorrelate on a timescale on the order of 3/Av hp . The 
MCF contains the peak from the autocorrelation of the 
intensity minus 2 /3 the square of the peak of the autocor- 
relation of the signal itself. This last contribution works as 
a partial bandpass correction, but not as a full bandpass 
correction, as the two contributions are not necessarily 
equal. 

3.1. Correcting for the bandpass 

The noise power is dominated by the sky background ra- 
diation at the frequencies used in this work. We therefore 
take the shape of the power spectrum of the signal in 
the off-pulse as an estimate of the bandpass. The WSRT 
bandpass contains modulation with a width on the order 



To make a clear distinction between time and frequency, 
we use time (t) instead of phase (cj>). 
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of 2 MHz and interference peaks (see also iKouwenhovenl 
120001 chapter 4)). Fig. shows the bandpass of on-pulse 
and off-pulse emission for one polarisation state of pulsars 
B0950+08 and B0329+54, averaged over 2,000 pulses. 

To correct for the bandpass, we calculate and apply 
a bandpass correction as follows. First, we estimate the 
average bandpass as the power of the off-pulse signal as 
a function of frequency averaged over all pulses. Wc then 
normalise by dividing each frequency value by the average 
of the bandpass. Before the calculation of the MCF, we 
Fourier transform the signal inside each window into the 
frequency domain, divide the value at each frequency by 
the square root of the corresponding value of the band- 
pass and Fourier transform back into the time domain. 
This process does not correct for a possible frequency de- 
pendent phase shift introduced by the system. We carry 
out the procedure described above separately for both po- 
larisations. High-level interference peaks are replaced in 
the frequency domain by the average of their left and 
right neighbours. Adjacent peaks are considered as one. 
Remaining interference peaks should only have a small 
influence on the MCF and only on long timescales. 





Fig. 2. Our results of the MCF calculated for PSR 
B0329+54 after bandpass correction. See Fig. 2] for ex- 
planation. 

In Fig0we see the corrected MCF for PSR B0329+54. 
The peak at low A<j> has gone down slightly, but remains 
present. Looking at Fig. |2t>, we see that there is strong- 
modulation in the on-pulse bandpass of PSR B0329+54 
which is absent in the off-pulse bandpass. This reduces 
the effectiveness of the bandpass correction, as this mod- 
ulation makes the on-pulse bandpass different from the 
off-pulse bandpass. We did not do a bandpass correction 
for PSR B0950+08, as there is no sign of coherence present 
in our data from this pulsar. 

3.2. Effect of scintillation on the MCF 

Having established that variations in the bandpass shape 
can theoretically influence the MCF, we now consider 
the additional influence of interstellar scintillation (ISS). 
Diffractive ISS introduces structure in the instantaneous 



spectrum of a pulsar and can therefore, in principle, influ- 
ence the MCF. In particular, we argue that ISS induces a 
signal in the MCF that is likely to masquerade as a false 
detection of temporally coherent pulsar radiation. 

Diffractive interstellar scintillation is observed in pul- 
sars, including B0950+08 and B0329+54, at centimetre 
wavelengths and lower. It imposes large (> 100%) tempo- 
ral and spectral modulations in the intensity of the radi- 
ation. The decorrelation timescale and bandwidth of the 
fluctuations are chiefly determined by two parameters: the 
diffractive scale rdiff, and the Fresnel scale rp. It is often 
convenient to consider the phase fluctuations imposed by 
the scattering medium confined to a thin phase screen of 
distance L from the observer. Then the diffractive scale is 
the distance on the phase screen over which the root mean 
square phase difference is one radian, and the Fresnel scale 
is set by the distance to the screen and the wavelength 
r F = y/XL/2-rr. 

Diffractive scintillation imposes random spectral vari- 
ations of characteristic bandwidth 



A^dc 



( r J¥L 



(7) 



across the pulsar signal, and this spectrum changes on a 
time scale ~ r^is/v, where v is the velocity of the scintilla- 
tion pattern across the Earth. The expected decorrelation 
bandwidth for observations of B0329+ 54 and B0950 +08 
at their o bserving frequency are 29 kHz <|Cordesll98rl and 
0.22 GHz <|Philli P s & Clegdll992|) . respectively. These val- 
ues were scaled from a frequency of 1 GHz and 51 MH z, 
respectively, assuming Ais^c oc v^ A l|Cordes et al.lll985h . 

Let us consider the effect of diffractive spectral varia- 
tions on the quantity Ci(At). The instantaneous power 
spectrum of the pulsar signal contains large variations 
with a characteristic bandwidth Av^c- Thus, since the 
amplitudes of the Fourier transformed voltages contain 
ripples, the observed voltages are broadened. In the case 
of scintillation, the broadening c an be described b y the 
pulse broadening function (PBF) j Williamsonll 1 972f) . The 
instantanous PBF is expected to have wiggles on a time 
scale set by the reciprocal of the decorrelation bandwidth, 
At « 1/A^dc- This implies that the intensity also con- 
tains variations on the same characteristic time scale. 
Neglecting any intrinsic temporal coherence due to the 
pulsar, the intensity autocorrelation function C/(Ai) is 
expected to peak at At = and decorrelate on a timescale 
set by the inverse of the scintillation decorrelation band- 
width. The average intensity autocorrelation function, ob- 
tained by combining the autocorrelations from many indi- 
vidual pulses and even over many diffractive timescales, is 
expected to exhibit the same decorrelation timescale. This 
is because every set of data is expected to exhibit spectral 
structure with a similar decorrelation bandwidth. 

The second contribution to the MCF comes from the 
autocorrelation of the signal voltages, Cs(At). The au- 
tocorrelation function of the pulsar voltages, S(t), and 
the on-pulse bandpass, S^Sfa)* are related by a Fourier 
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transform: 



Cs(At) ^ (S(v)S*(v)), 



(8) 



where the angular brackets denote an average over many 
scintles in either/both the temporal and spectral domains. 
Thus the autocorrelation Cs(Ai) is the inverse Fourier 
transform of the mean bandpass. However, the mean band- 
pass contains little structure due to the stochastic scin- 
tillation. Observations over many scintillation time scales 
wash out the frequency dependent structure due to scintil- 
lation, so that the mean on-pulse bandpass reflects mainly 
the instrumental response to the intrinsic spectrum of the 
pulsar. Thus scintillation has little effect on Cs(Ai). 

3.3. The Modified Coherence Function due to 
Scintillation 

A set of synthetic observations from data generated by 
a scattering simulation is constructed in order to demon- 
strate the effect of scintillation on the MCF suggested by 
the preceeding arguments. 

The simulations are conducted by constructing a set 
of phase fluctuations as a function of position and wave- 
length, </>(r,A), frozen onto a thin phase-changing screen 
located a distance L from the observer's plane. We take a 
plane wave of unit amplitude incident on the phase screen, 
so that the phase of the wave upon exiting the screen is 
4>{r). The phase fluctuations are generated according to a 
von Karman power spectrum l|Goodman fc Naravanll 



q) = Qo [ |q 



-0/2 



exp(-|q|/<?max), 



(9) 



where Qq is the amplitude of the power spectrum and 
sets the diffractive scale length, rdiff. The Kolmogorov 
valu e of = 11/3 is assumed for the power law in- 
dex I Lee fc JokiDiill975h , and g m ax and q min are the upper 
and lower cutoffs of the power spectrum respectively. Each 
specific realisation of the scattering screen is generated by 
randomly choos ing the phase associated wi th each spec- 
tral mode (e.g. iGoodman fc Naravanl |l989), so that the 
Fourier transform of the phase fluctuations, </>(x), is given 

by 



q) = Q0(q) 1/2 ex P [i|q|-x] 



(10) 



with X £ (0,27r]. The wavefield at the observer's plane 
is the inverse Fourier transform of the Fourier transform 
of the wavefield at the exit plane of the phase screen, 
exp[i0(r)] multiplied by the Fourier transform of the prop- 
agator exp[— i{2irr-p) 2 q 2 /2\. 

Temporal fluctuations in the wavefield are obtained by 
moving the phase screen relative to the observer, however, 
due to the limited number of grid points short time scale 
variations could not be simulated. Spectral variations are 
obtained by scaling the phase fluctuations according to 



M) 



A 



(r,A ), 



(11) 



and by altering the Frcsncl scale in a similar manner, 
rp(A) = \A\AAo r F(Ao) where A is a wavelength within the 
observation bandwidth and Ao is an arbitrary, but fixed 
wavelength. Synthetic pulsar observations are generated 
by shifting the wavefield due to the scintillation pattern 
in time and frequency in this manner, and by choosing 
the length scales rdiff and r-p appropriate to the scatter- 
ing conditions. The MCF is computed directly from the 
wavefield, where the wavenumbers, fc m ; n and fc max , were 
chosen to match a 5 MHz bandwidth around 328 MHz and 
the upper and lower cut-offs of the power spectrum, <7 m i n 
and q m ax, were set at values corresponding to a physical 
size much larger than the scale of the phase screen and 
much smaller than the size of one discrete point, respec- 
tively. Table El shows the values of the other scintillation 
parameters with a brief explanation. The distance of the 
scintillation screen, D, was chosen to be in the middle be- 
tween the pulsar and the observer. By tuning Qo, we tried 
to create scintillation similar to that in our data from PSR 
B0329+54. The simulation results for PSR B0329+54 are 
shown in Fig. 0]and clearly show an excess similar to that 
in our MCF of PSR B0329+54. This simulation does not 
reproduce the peak in the first bins. However, in order 
to see the effects of the short time scale variations of the 
scintillation pattern, we have simulated a growing scintil- 
lation pattern similar to that visible in Fig. Eh> on top of 
a flat bandpass. Note that the individual modulations in 
the bandpass of our data from PSR B0329+54 often over- 
lap each other. This simulation did result in a peak in the 
first bins of the MCF. Furthermore, we split up our data 
from PSR B0329+54 into smaller sets of data, calculated 
the MCF of each dataset and averaged them together. We 
did this for datasets ten to hundred times shorter than 
the total length of our data and saw that the peak in the 
first bins of the MCF decreases when the dataset becomes 
smaller. 

4. Discussion and Conclusions 

We find that frequency modulation due to scintillation and 
possibly the shape of the WSRT bandpass have an effect 
on the MCF at small delay values. This may be problem- 
atic in determining a coherence time. For PSR B0950+08 
we find no signature of coherence. For PSR B0329+54 we 
find two features at small delay values: a peak in the first 
two bins and a broad excess reaching up to 35 /xs. We have 
shown that the shape of the bandpass can theoretically 
cause a peak in the first bins of the MCF. The bandpass 
correction, which uses the off-pulse, might not be effective 
in the case of PSR B0329+54 due to the frequency mod- 
ulation present in the on-pulse. By comparing Fig. 0] with 
Fig-E 3 , we see that scintillation reproduces the broad rise 
observed in the MCF of PSR B0329+54. Furthermore, a 
simulation of an increasing scintillation pattern on top of 
a flat bandpass, where the modulation due to the scintil- 
lation was overlapping, has shown that small changes in 
the scintillation pattern can cause a peak in the first bins 
of the MCF. This is similar to a time dependent variation 
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Fig. 3. Bandpass of one polarisation of WSRT for (a) On-pulse of B0950+08, (b) On-pulsc of B0329+54, (c) Off-pulse 
of B0950+08, (d) Off-pulse of B0329+54, all averaged over 2,000 pulses. The peaks are caused by interference. The 
on-pulse of B0329+54 shows scintillation with a width of approximately 30kHz (or 25 bins). 



in the bandpass itself. Moreover, when smaller timescales 
are used to calculate the MCF the peak is seen to decrease 
also indicating that scintillation may play a role here. 

From the above, we conclude that scintillation can be 
responsible for both the broad excess as well as the peak 
in the first two bins of the MCF of PSR B0329+54. 

We now discuss whether the excess found in the MCF 
of three pulsars by JAP could also be due to scintillation. 
In Table E3 we show the similarity between the diffrac- 
tive decorrelation bandwidth of the different pulsars and 



their coherence time, defined as the point where the MCF 
becomes zero. The values for Av^c were scaled from a fre- 
quency of 51 MHz for PSR B0950+08 and fr om I GHz for 
the o ther pulsars, assuming Av ^ oc v iA llCordes et ID 
Il985h . No errors were quoted in lCordesl <)l98fil) . however 
it is known that Avd c can vary signific antly with time, i n 
some cases as much as a factor of 2—3 l|Bhat et al.l ll999). 
Assuming that the rise in the MCF is due to scintilla- 
tion, we estimate the MCF to become zero when A(f> is 
on the order of l/Av^c- Looking at Table we see that 
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Table 3. Diffractive decorrelation bandwidth and coherence time of pulsars B0329+54, B0823+26, B0950+08, 
B1133+16 and B1937+21. 



Pulsar 


Frequency 
(MHz) 


Time resolution 
(ns) 


A^dc 


1/A^ dc 
(Ms) 


Measured coherence 
time (/is) 


B0329+54 


328 


100 


29 kHz [1] 


34 


35 [4] 


B0823+26 


430 


100 


0.81MHz [1] 


1.2 


1.5 [5] 


B0950+08 


382 


50 


0.22 ± 0.03 GHz [3] 


4.5xl0~ 3 


none [4] 


B0950+08 


430 


100 


0.35 ± 0.05 GHz [3] 


2.9xl0~ 3 


0.4 [5] 


B1133+16 


430 


100 


1.47 MHz [1] 


0.67 


1.1 [5] 


B1937+21 


430 


100 


4.2 ± 0.9 kHz [2] 


238 


none [5] 



[1] Cordes (1986) 

[2] Cordes et al. ( 19901 

[3] Phillips fc Clegg Q992) 

[4] Our data 

[5] Data from JAP 



Table 2. List of scintillation parameters and values used 
in the simulation 





Input 

parameters 



Discription 



Fig. 4. The MCF of simulated scintillation at a frequency 
of 328 MHz and a bandwidth of 5 MHz. The input param- 
eters are set to generate scintillation similar to that in our 
data from B0329+54 



the values of the reciprocal of the diffractive decorrela- 
tion bandwidth and the measured coherence time (fourth 
and fifth column, respectively) are indeed similar. For our 
result of PSR B0950+08 there is no measured coherence 
time, while JAP find a coherence time for this pulsar of 
0.4 u s. Assuming the diffract ive decorrelation bandwidth 
from IPhillios fc Cleggl lll992l) the time resolution of our 
data and that of JAP are too large to see the scintillation. 
The observed feature in the MCF for PSR B0950+08 of 
JAP (see their Fig. 1) might then result from a time varia- 
tion in the bandpass of their data. There is, however, some 
controversy as to what the diffractive deco rrelation band- 
width for PSR B0950+08 is. According to lCordesI l|l986|) 
A^jc is 4.0 MHz at an observing frequency of 430 MHz. 



fcmin = 6.807 m 1 
fcmax = 6.912 m -1 



D = 0.7 kpc 

Q =l.lxl0- 41 m" 17 / 3 
N = 128 



Minimum wavenumber 
Maximum wavenumber 
(fcmin and fc max are the edges of 
the observed frequency range) 
Distance between scintillation 
screen and observer 
Amplitude of the power spectrum 
of the phase fluctuations 
Number of discrete points across 
the wavefield (Total number of 
grid points — N 2 ) 



Measured 
parameters 



r F = 


1.7 x 10 9 m 


Fresnel scale (= -JD/k) 


rdifi = 


= 3 x 10 7 m 


Diffraction length scale 


J"rcf = 


= 1.7 x 10 11 m 


Refraction length scale 






(= r F 2 /r dift ) 


Ai/dc 


= 0.02 MHz 


Diffractive decorrelation 






bandwidth 



This gives a value for l/Ai/dc of 0.25 /is, which is in the 
order of size of the coherence time of 0.4 /is, measu red 
by JAP. Assuming the value given bv ICordea l|l986|) . it 
would appear that the signal-to-noise ratio in our data is 
insufficient to see the scintillation, as we see no excess in 
the MCF. For PSR B1937+21, JAP also found no coher- 
ence time. Although this pulsar is known to scintillate, the 
value for l/Ai^c is so large that we can expect the excess 
in the MCF due to scintillation to be smeared out over 
such a large range of delays as to make it unmeasurable 
with the sensitivity of their observation. 

We conclude that scintillation and possibly the shape 
of the bandpass causes the excess in the MCF of our data 
from PSR B0329+54. Furthermore, we conclude that the 
coherent features, found by JAP, also appear to be the 
result of scintillation. We therefore cannot confirm that 



8 



J.M. Smits et al.: On the Search for Coherent Radiation from Radio Pulsars 



the MCF is clearly showing us the presence of coherence 
in these pulsars. 

The authors would like to thank J. Cordes for his ex- 
tensive comments and F. A. Jenet for his helpful discus- 
sions which have both greatly contributed to the accuracy 
and clarity of this paper. 
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